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Abstract 

Classical-quantum correspondence for conservative chaotic Hamiltonians is investigated in terms of the structure of 
the eigenfunctions and the local density of states, using as a model a 2D rippled billiard in the regime of global 
chaos. The influence of the observed localized and sparsed states in the quantum-classical correspondence is 
discussed. 
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The subject of this Letter concerns the quantum-classical correspondence of autonomous Hamiltonian 
systems with classical chaotic dynamics. As a paradigm for this class of systems we consider the motion 
of particles in the 2D rippled billiard depicted in Fig. 1. It is known [1-3] that the dynamics of a classical 
particle in this billiard undergoes the generic transition to chaos (regular-mixed-global) as the amplitude 
a is increased (d is fixed). Thus the results that can be obtained by studying this particular model are 
applicable to a wide class of systems. A modern physical realization of this rippled billiard of finite length 
is a mesoscopic electron wave guide. In [1] its classical transport properties, such as resistivity, were related 
to its dynamical properties, yielding a transport signature of chaos. The analysis of the quantum motion in 
the infinitely long rippled channel (the periodic rippled channel) is useful for the understanding of universal 
features of electronic band structures of real crystals [4] , propagation in periodic waveguides [5] , quantum 
wires [6-8] and films [9]. Signatures of chaotic diffusion in the band spectra of a similar periodic billiard 
have been investigated in Ref. [10] . In [2] the energy band structure of the quantum version of the periodic 
rippled channel was calculated and certain aspects of the quantum-classical correspondence were investigated 
in terms of the Husimi distributions. 



y = d + acos(x) 








FIG. 1. Geometry of the rippled channel. 



The standard signature of chaos in quantum mechanics [11] based on the random matrix conjecture [12] 
showed a qualitative agreement for the energy level spacing statistics of the periodic rippled channel [13]. 
However, a more detailed study [14] revealed important deviations. Motivated by this discrepancy, which 
may be found in other chaotic systems, here we use a novel approach that has been suggested in Ref. [15] 
to analyze the quantum-classical correspondence of complex systems. Such an approach relies on the study 
of the structure of the local density of states (LDOS), and the shape of the eigenfunctions (SEF) in the 
basis of non-interacting particles (quasi-particles) [16]. It has been confirmed that the LDOS and SEF 
have well defined classical counterparts in the chaotic regime in dynamical systems such as that of two 
interacting spin particles [17] and the three orbital schematic shell model [18], see also very recent works in 
[19]. Below we extend these studies to chaotic billiards by treating the degrees of freedom as independent 
particles and incorporating the effects of the boundaries into the Hamiltonian operator. We show how to 



construct the classical counterparts of LDOS and SEF for chaotic billiards. We shall demonstrate that these 
quantities determine the global shape of the quantum LDOS and SEF, and that the phenomenon of quantum 
localization within the energy shell manifests itself as strong fluctuations about the global shape. We remark 
here that this localization occurs not in configuration space but in energy (or momentum) space (for chaotic 
billiards this kind of localization was first investigated in Rcf. [20]). 

We limit our studies here to the case of global classical chaos which occurs in wide [d > 27r) channels [21]. 
The control parameter that determines the degree of chaoticity, obtained by the standard linearization of its 
mapping, is given by K = 2 -^ 2; . In agreement with Chirikov's overlapping criteria our numerical experiments 

show global chaos for K ~ 1. We have compared the phase space dynamics for various combinations of a 
and d yielding the same value oi K > 1, and found their Poincare maps to be indistinguishable from each 
other. In what follows we take (a/2n, d/2n) = (.06, 1), i.e., K w 1.5. 

The quantum mechanical description involves the solution of the Schrodinger equation H = \{p 2 z + Py) 
subject to the boundary conditions ip(x,y) — at y = and y = d + acosx. It is both convenient and 
illuminating to go to the new coordinates u — x, and v = d+ ^, s (x) ' ^ or * n these coordinates the boundary 
conditions become simpler: *&(u, v) — at v — 0, d. On the other hand, the Hamiltonian acquires a much 
more complicated form, 

i?=f^ 1/4 ^5 Q V /2 ^. 9 - 1/4 , a,P = u,v, (1) 

which is simply the kinetic energy expressed in covariant form [22]. The momentum is now given by P a = 
—ih[d a + jd a ln(g)] — —ihg~ 1 / 4 d a g 1 / 4 , where g af3 is the metric tensor, and g = Det(g a p) = [1 + ecosu] 2 , 
e = a/d, is the metric (see [2] for details). Thus, we have formally transformed the model of one-particle 
motion in the rippled billiard to that of two interacting particles (identified with the two degrees of freedom). 
In this way, the complexity of the boundary in the original model is incorporated into the interaction 
potential. 

Since the Hamiltonian is periodic in u, the energy eigenstates satisfy Bloch's theorem: ^e{u, v) = 
cxp(z/cu)<i>fc(u, v) with $fe(u + 2ir, v) = "^(u, v) and Bloch wave vector k = k(E). For an infinite channel, k 
takes a continuous range of values and we choose the first Brillouin zone to lie in the interval [— \ < k < i]. 
Further, we can expand v) in a double Fourier series and write the a th eigenstate of energy E a (k) as 

oo oo 

* a (u,v;k)=J2 E CZ n (k)ct> k mn (u,v) (2) 

m—1 n— — oo 

where 

4„(«, v) =<u,v\m,n > fe = 7r- 1 / 2 ff - 1 / 4 sm(m7r«)e J ( fe+n ^. (3) 

The factor n^^g- 1 / 4 comes from the orthonormality condition in the curvilinear coordinates (u,v). For 
future reference we shall call m the mode (or channel) number and n the Brillouin Zone number. Note 
that <j> mn are the eigenstates of the unperturbed momenta squared (and, therefore, of the unperturbed 
Hamiltonian): P^(j)^ nn = h 2 '(rm:) 2 (j) mn , P 2 <t> mn — h 2 (k + n) 2 (j)^ nn . The problem requires the diagonalization 
of the matrix Hiy (k) =< I \ H \ I' where | I >k=\ m, n Note that for each pair (m, n) we associate a 
number /. Even though the energy spectra is independent of the choice of the function l(m, n), the structure 
of the matrix clearly depends on it. For our purposes, it is essential (as it will become clear below) that the 
unperturbed Hamiltonian basis be ordered in increasing energy: E® +1 {k) > E^{k), where H°(j)f — E[(k)<f)f, 

Ei(k) = \((n+k) 2 + (m7r) 2 ). This defines the choice of assignment (m, n) => I. Once the matrix Hiji(k) has 
been diagonalized, its (exact) eigenstates ^ a (k) = ^2 C"(fc)0; are also re-ordered in energy (E a+1 > E a ). 
We adopt the convention that the Greek supcrindcx (Latin subindex) denotes the exact (unperturbed) state. 
The amplitudes C"(fc) form the state vector matrix; its elements along the row a are the components of 
the a th state in the representation of the unperturbed re-ordered basis, and the elements along the column 
I give the unperturbed state / expanded in the re-ordered perturbed basis. We shall study the structure of 
both, rows and columns, the latter being essential in the construction of the local density of states (LDOS) 
of the system. 
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In the previous work [13] wc studied the level spacing statistics ol the rippled channel in connection with 
the Random Matrix Theory conjecture. Using the Bloch momentum k as an external parameter we found 
that the level spacing statistics is the same for all values of k (except k « 0) within the Brillouin zone. Thus, 
here we arbitrarily take the value k = 0.1 to explore the structure of the exact eigenstates and from now on 
we drop, for economy of notation, the label k. 
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FIG. 2. Typical high energy eigenfunctions for a) a = 2078, b) a = 2079 and c) a — 2080. The localization 
measures are: a) Ih = 653.5, h pr = 411.5; b) Ih = 4.5, h pr = 3.6; and c) Ih = 747.8, h pr = 571.7. The amplitude 
and width of the channel are a/2n — .06 and d/2n — 1, respectively. The value of the Bloch wave vector is k = 0.1. 

Fig. 2 shows three consecutive high-energy eigenstates (C™, a=2078, 2079, and 2080) in the re-ordered 
unperturbed basis. This sequence illustrates the kind of eigenstates that can occur typically; in a given 
sequence for this range of values of a roughly 15 percent are states such as the one shown in Fig. 2b, and 
about 60 percent are like those of Figs. 2a and 2c. The rest are intermediate between these two types: 
"sparse" states characterized by a few dominant peaks in a background of smaller components. Eigenstates 
can be characterized by various localization measures, such as the "entropy localization length" l-u and the 
"inverse participation length" l ipr (see, e.g., Ref. [17]) defined by 

l n = 2.08cxp(-H), l ipr = 3/P, (4) 

where 

N N 

H = J2 I r(E) T In | r(E) | 2 , and P = J2\ * a {E) | 4 . (5) 

a—l a—1 

These two measures are found to be two orders of magnitude smaller for the eigenstate 2079 than for the 
eigenstates 2078 and 2080 shown in Fig. 2. Extensive numerical results [23] show that these lengths fluctuate 
wildly as a function of the number of the eigenstate a, or, equivalently, as a function of its corresponding 
energy E a . As for extremely localized states such as that of Fig. 2b, they are expected to play an important 
role in transport since they are practically the same as the unpertubed (regular) states of the model. In the 
energy spectrum, they are located at the edges of the cells defined by different values of the Brillouin Zone 
number n. Moreover, these states have the largest group velocity dE/dk; see [23] for details. 

Our main interest here is to relate the structure of the matrix C" to properties of the corresponding 
classical system. Since C" =< | <f>j > as a function of j (or equivalently, E°j) is the projection of 
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the exact state onto the states of the unperturbed system, a classical counterpart of | C" | 2 as a function 
of j can be obtained [15] by projecting the (chaotic) trajectories generated by the Hamiltonian H onto the 
unperturbed Hamiltonian H°, where H = H° + V and V is the perturbation. Specifically, this can be done by 
substituting the orbit <&(£) = (x(t),y(t),p x (t),p y (t)) of a typical initial condition generated by H with energy 
E onto H°. Then the energy of the unperturbed Hamiltonian varies in time, E°(t), fluctuating within a range 
AE which defines the width of the so-called "energy shell." After a sufficiently long time a distribution of 
energies W c i(E°/E) can then be obtained from E°(t) for a single (chaotic) orbit. Alternatively, W c i(E°/E) 
can be calculated for shorter times but taking several orbits. The equivalence between these two ways is 
expected if motion generated by the full Hamiltonian is globally chaotic. 

The above discussion outlines the general procedure that claims that W c i (E / E) is the classical counterpart 
of | 2 as a function of E®. It has been applied [15-19] in a straight forward manner in systems 

where H° and V can readily be identified. For chaotic billiards (with Neumann or Dirichlet boundary 
conditions) the Hamiltonian of both unperturbed and perturbed systems is simply the kinetic energy and 
the non-intcgrability comes from the boundary conditions. In order to incorporate the perturbation into the 
Hamiltonian operator, it is necessary to perform an appropriate coordinate transformation in the same way 
as it was done above for the quantum model. To accomplish this, we take the quantum Hamiltonian, Eq. 
(1), and commute all the operators, transforming them to c- numbers, 

H = \g aP P a P P , (6) 



subject to the "flat" boundary conditions in (u, v) coordinates. Expanding this expression, we obtain 

(7) 



" l + ei u v+ (1 + eO 2 



where £ = cosu, £ u = d^/du, and e = a/d. The canonical transformation from the old to new variables is 
given by u = x, v = yd/(d + a£), P u = P x + ey£ x P y /(l + e£), and P v = (1 + eQPy We now separate the 
Hamiltonian H = H° + V analogously to the quantum model: H° = \ [P 2 + P 2 ] ; V = H - H°. 

To compare the classical distribution W c i(E° / E) with the quantum one | C" | 2 (as a function of E°j) we 
take the energy of the classical Hamiltonian to be equal to the energy of the state | \& Q > under consideration. 
Since the classical phase space dynamics is independent of the energy of the particle (the kinetic energy can 
be rescaled) we arbitrarily take it to be equal to one and introduce an effective Planck's constant U e f / . 
In Fig. 3 we show the classical and quantum distributions (shapes) of four eigenstates. The areas under 
the classical and quantum curves have been normalized to 1. These states were chosen from four different 
energy regions (different values of Ti e ff) an d they represent typical states in each region of about one hundred 
eigenstates. These are typical in the sense that most states in each interval are characterized by typical 
localization lengths in that range. For example in the region of 700 < a < 800, the length In oscillates 
around In = 200 with only 8 eigenstates with In < 50 and the maximum of In is about 340. That is, the 
states of Fig. 3 are neither extremely delocalizcd nor extremely localized. 

A simple inspection of Fig. 3 indicates that for low energies few unperturbed components participate in 
the construction of the exact eigenstate and as the energy increases, the participation of more and more 
components increases. Based on the Weyl formula and the condition that the de Broglie wavelength be equal 
to the amplitude of the ripple, one obtains an estimate for the number of levels defining the perturbative 
border N p . For the geometrical parameters of the channel (a/2n = .06, d/2ir = 1) we have N p w 139. 
Thus, for low energies, as exemplified by Fig. 3a, there is hardly any resemblance between the classical and 
quantum distributions, since at such energies the system is still in the perturbative regime and far from 
the semiclassical domain. It is also clear that as the energy increases the global shape of the quantum 
distribution is closer to the classical distribution. It is remarkable that the range of energy components that 
participate in the construction of the perturbed eigenstate can be deduced from the classical distribution. 
This range of energies defines the energy shell of the eigenstate under consideration. Note that the slight 
assymmetry observed in the range of the quantum distribution with respect to E° = 1 is the same as in the 
classical distribution and can be explained analytically, see details in Rcf. [23]. 
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FIG. 3. Individual eigenfunctions | C" | 2 in the energy representation and the classical distribution W c i(E° /E) 
(thick curve) as a function of the energy E°j for a— a) 101; b) 736; c)1361; and c)2000. The entropy localization 
length l n is a) 18.7; b) 162.7; c) 422.3; and d) 707.2. The effective Planck's constant h eff is a) 0.24, b) 0.09, c) 
0.07, and d) 0.05. Same values of amplitude a, width d, and k as in Fig. 2. Note the different vertical scales due to 
normalization conditions. 



For states with largest values of U pr or 1-h., our numerical experiments show a much better quantum- 
classical correspondence. This occurs because as more components are needed to form the perturbed state, 
then the number of components with amplitudes above the classical shape decreases. Even so, strong fluc- 
tuations of the quantum distribution with respect to the classical shape do not dissappear completely even 
for the case of extremely delocalized (within the energy shell) states, see Fig. 3d. 

We now try to smooth these fluctuations by calculating the averaged distribution, also known as the 
average shape of the eigenfunction (SEF) 

W(E 3 /E a ) = I C°j I 2 6{E - E a '), (8) 

a' 

where the sum is taken over a number of eigenstates in the neighborhood of the exact state >. 
Another quantum quantity that shall be compared with its classical counterpart is the so called "local 
density of states" (LDOS) or "strength function" , widely used in solid state and nuclear physics. It is defined 

as 

u(EyE° j ) = J2\C a j ,\ 2 5(E°-E j ,), (9) 

3' 

where the summation is done over a number of perturbed states in the neighborhood of the unperturbed 
state | j >. This quantity gives information about the evolution of the system. In particular, the spreading 
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width of this function determines the energy range associated with its "decay" into other states due to 
the interaction. The classical counterpart of the LDOS is constructed in much the same way as for the 
"classical eigenfunction" distribution W c i(E°/E) except now we project the trajectories of the unperturbed 
Hamiltonian H° onto the exact Hamiltonian H [24] . 
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FIG. 4. a) Shape of eigenfuctions (SEF) in the energy representation W(E°j/E a ) and its classical counterpart 
W c i(E°/E) (thick curve); b) LDOS ui{E a / E°j) and its classical counterpart u> cl (E/E°) (thick curve). Same values 
of amplitude a, width d, and k as in Fig. 2. For the SEF the average is taken over the interval 1870 < a < 1950 and 
for the LDOS the average is over the interval 1900 < j < 1990. 



Figures 4a and 4b show both distributions, the average SEF and the LDOS, together with their classical 
counterparts. Even though the agreement between classical and averaged quantum eigenfunction shapes has 
improved compared to the individual shapes of Fig. 3, the quantum fluctuations are still relatively large. 
Moreover, the averaging procedure reveals a three-peak structure for both, classical and quantum models. 
The origin of these peaks can be undestood by a detailed analysis of the classical trajectories. In particular 
for SEF, it was found [23] that the chaotic trajectories dweel for a long time in the neighborhood of stable 
and unstable period-one orbits. The unstable (stable) period-one orbit is defined by P u — and x — 0(ir). 
The right peak corresponds to motion perpendicular to the x axis at x — tt (unstable orbit). Similarily, 
the left one results from the stable orbit at x = 0, see details in Ref. [23]. In contrast, the central peak is 
produced when the trajectory is parallel to the x-axis. A similar analysis explains the origin of the structure 
of the classical LDOS [23]. One can detect a difference in the structure of the LDOS in comparison to 
that of the SEF, mainly, in the central region. This fact may be explained by the difference between time 
averages and phase space averages for finite times (see also discussion in Ref. [25]). It is remarkable that the 
quantum distributions are sensitive to the existence of these periodic orbits. In general there is an overall 
good correspondence between classical and quantum distributions (as exemplified by Fig. 4) but with large 
fluctuations due to the existence of localized and sparse states within the energy shell. 

The kind of localization which we observed in this study occurs in a finite energy range determined by 
the classical energy shell. On the other hand, in time-periodic models, localization happens in unbounded 
momentum or energy space, and chaotic eigenstates are dense in the unperturbed basis. On the contrary, 
localization in conservative systems reveals itself in two different ways, one is very strong (perturbative) 
localization, see example in Fig. 2b, and the second one appears in the form of sparse eigenstates, see Fig. 
3a. So far, it is not clear whether in conservative systems, the former type of localization can also appear. 

In summary, we have analyzed quantum-classical correspondence for the eigenstates and the local spectral 
density of states for non-integrable billiards, using as a paradigm a 2D channel with periodic sinusoideal 
boundaries. We have shown that the global shape of the quantum quantities can be determined from their 
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classical ones. In particular, the range of the energy shell and its global shape can be deduced from the 
classical distributions. It was also shown that the existence of localized and sparse states within the energy 
shell give rise to strong quantum fluctuations. These states are also responsible for the observed deviations 
from random matrix theory predictions for level spacing statistics and are expected to play an important 
role in the transport properties of the system. We remark that the comparison was done here for quantum 
and purely classical quantities. Thus, it is very interesting to develop a kind of semiclassical approach for the 
shape of eigenstates and the LDOS, based on the trace formula (for time-periodic pertubation, this problem 
has been addressed in Ref. [26]). 
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